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abstract 

Energies and wave functions of edge states in two dimensional electron 
evaluated for a finite step potential barrier model. The spectrum, 
instead of smooth bending of Landau branches in the vicinity of the barrier 
acquires a steplike form; unexpected edge plateaus and significant energy 
gap reduction take place between the neighbouring Landau branches above 
the barrier tops. The origin of these phenomena is traced. Stability with 
respect to modifications is established. Manifestation of the qualitatively 
new features of electron densities of states in abrupt confinements through 
magnetooptical and nuclear magnetic spin relaxation effects is proposed. 
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1 Introduction 



Strong magnetic fields and progress in heterostructure quality allowed dis- 
covery of such dramatic phenomena in two dimensional electron gas (2DEG) 
like Integer and Fractional Quantum Hall Effect [?]. In the microscopic the- 
ory the edge states, classically represented by electrons skipping in circular 
segments along the edges, play an important role. Many magnetotransport 
experiments in two dimensional electron gas have been qualitatively under- 
stood recently by means of a simple edge - state model [?]. This model is 
based on the picture of smooth Landau level bending by the potential formed 
by external charges. The intersections of each Landau level with the Fermi 
surface create widely separated narrow edge channels [?]. 

The semiclassical notion of the Landau level bending, however, is not 
always applicable to the electronic structure in the vicinity of the barrier 
representing an edge or an interface (or random potential). It is definitely 
not valid for potentials with large gradients. In such situations often used 
quantum mechanical model [?] is limited to the extreme case of an infinite 
barrier. Recently Chklovskii et al [?], using the self-consistent electrostatic 
approach, showed that the resulting effective potential should acquire a step- 
like shape even if the external potential is smooth. 

In this paper, finite step potential is considered to better approximate 
the interface or the steplike potential just mentioned. The electron structure 
exhibits two rather unexpected phenomena: the edge plateaus and the Lan- 
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dau gap reduction. These phenomena are not present in either of previously 
studied models. 

The interface might be either a "boundary" confining 2DEG or an inter- 
face between two different materials. The finiteness of the barrier confining 
2DEG is especially important, if the confinement is realized by means of an 
interface between two similar materials. Then, in very strong magnetic fields, 
the magnetic field induced splitting (between Landau levels) may be compa- 
rable to the potential barrier height. For example, in the 2DEG formed at 
the interface between GaAs and GaAl 1 _ ;E As^, the effective interface poten- 
tial barrier is about 0.3 eV. The Landau level spacing in strongest magnetic 
fields experimentally available is just a few times smaller. In particular, for 
the steplike potential with wide steps (on the scale of magnetic length) of 
height hu c [?, ?], the two quantities are equal. 

Two dimensional electron systems are realized at abrupt interfaces be- 
tween two slightly different semiconductor crystals. Such systems are pre- 
pared by molecular beam epitaxy where deposition as well as doping can 
be controlled on atomic layer level. Lateral confinement within the plane of 
2DEG can be achieved in various ways. In one of them, the gate voltage 
restricts the electron motion. Due to the distances between the gate and the 
2DEG plane, the resulting confinement potential is very smooth within the 
plane. Mesa etched samples in which the surface charge substitutes the gate 
voltage, fall also into this category of smoothly confined systems. These have 
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been studied in detail theoretically [?]. 

Recently, systems with varying chemical composition along the growth 
plane have been prepared. In such systems the effective confinement barrier 
height is approximately given by the well studied discontinuity in valence and 
conduction bands at heterostructures. Moreover, these barriers are localized 
within just a few interatomic distances. Therefore on a scale of magnetic 
length, even for strong magnetic fields they represent very abrupt barriers. 
For such systems our model is appropriate. Systems of this type have been 
also prepared in a form of very narrow channels (quantum wires). In these 
systems consequences of interaction between the two opposite edges can be 
studied. 

The paper is organized as follows. In section 2, general properties of con- 
finement barriers are reviewed. In particular, the concept of lagging of the 
electron center of mass behind the Larmor orbital center, due to the edge, 
is introduced. Next section is devoted to a rather detailed analysis of the 
simplest model of an abrupt barrier: the finite rectangular step. Spectrum 
exhibits two rather unexpected phenomena. Landau branches do not rise 
steadily towards the edge, but instead develop a series of "edge plateaus". 
The energy separation between the Landau branches does not remain con- 
stant: at certain regions within the edge it gets significantly reduced com- 
pared to Tiuj c . Both of these phenomena are peculiar to finite abrupt edge 
barriers and do not exist in smoothly confined systems. In section 4, some ex- 
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periments in which the characteristically modified electron densities of states 
should be observable like nuclear magnetic spin resonance and magnetoop- 
tics are discussed. The stability of results obtained for the simple model is 
examined in section 5. 

2 General description of the edge. 

Here, we consider a system of 2D noninteracting electrons in the vicinity of 
a boundary under the homogeneous magnetic field B perpendicular to the 
xy - plane. The edge or interface is described by a potential barrier V(x). 
The one particle Hamiltonian is 

where A(x) is the vector potential describing the applied magnetic field. In 

— * 

the Landau gauge, A = (0, Bx, 0), the motion along the y - direction is free 
and we can separate variables: 



i>n,x(x,y) = 




The wave function 4>n,x{x) obeys the one dimensional Schrodinger equation: 

<j> n , x (x) = E n {X)(t> n , x {x) (3) 

where I = \J ^ is the magnetic length, oj c = ^ is the cyclotron frequency 
and X is coordinate of the center of a Larmor orbit. This is the Schrodinger 



K 1 d 2 1 

- — -rr; + o^cOs - x ? + V(x) 
2m dx 2 2 
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equation for the harmonic oscillator superimposed with the barrier. The 
integer n parametrizes discrete Landau levels, n — 0,1,2, .... 

If the barrier V(x) is a smooth function of x (so that the force dv ^ is 
small compared to fiuj c /l), then all the Landau levels follow the potential: 

E n (X)= (n+^)hu c + V(X) (4) 

In this case the spacing between different Landau branches always remains 
a multiple of tuv c . 

In another extreme case of an infinite step barrier, simple boundary con- 
dition of vanishing of the wave function is to be imposed: <f> n ,x(0) = 0. 
It was studied by McDonald and Streda [?], who obtained energies of first 
few levels as a function of the distance X. Deep inside the region to the 
left of the barrier (X << 0) the influence of the interface is negligible and 

As X approaches the barrier, the energy levels rise 
due to repulsive effect of the infinite barrier. For orbits centered in the 
"forbidden" region to the right of the barrier, the energies continue to rise 
indefinitely. 

Certain general statements about the electron structure can be made 
without specifying the shape of the confinement barrier V(x). The Hamilto- 
nian of the ID Schrodinger equation depends on X. The Hellmann - Feynman 
theorem [?] applied to X as a parameter provides the following useful relation 
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which is independent of the barrier profile: 

dE n (X) ,dH x | , 2 i iv u /v\ 

=< <PnX\—jj^\<PnX >= mUJ c < <p n x\X ~ X\<p nX > (5) 

An important quantity X — x, where x is the average electron position, will 
be referred to as a displacement of the classical Larmor orbit center due to 
repulsion by the barrier. The displacement therefore becomes: 

mu>* dX 

The center of mass x characterizing the wave function <p n x{x) is thus easily 
obtained by differentiating the dispersion relation £J n (X)0. It is clear, that 
far from the barrier, x « X, the displacement is zero, and X = x. When 
the Larmor orbital center X approaches the barrier, its repulsion gives rise 
to lagging of x behind X. For the infinite confinement barrier, lagging of x 
behind its Larmor orbit's center X steadily increases as X grows (the wave 
function is nonzero for x < only). 

For any monotonously rising potential the function E n (X) is monotonous 
in Xf\ In the particular case of the finite step barrier, x laggs behind X (it 
will be proved later that for monotonously increasing V(x) > the deriva- 
tive > for any X and n). The magnitude of the displacement, 
however, for finite barriers, may rise and fall with growing X, as we will see 
in subsection 3.6.. 



This formula was derived by O. Heinonen and S. Taylor using variational principle, 
see [?]. 

3 This follows from a simple variational argument. 



3 The elementary rectangular step model. 

Let us now concentrate on the case of abrupt potentials. The basic potential 
of this kind is a single rectangular barrier V(x) = V6(x). Generalizations to 
several steps or to less abrupt barriers will be considered later. We start with 
an analytic perturbative treatment of small barriers which will be shown to 
contain all the qualitative features of the exact solutions. 

3.1 Approximate calculation. 

If V « Tiuj c or if a state is not very close to the interface, we can use the 
perturbation theory around Landau levels to calculate the dispersion relation 
E n (X). The perturbed energies to the first order in the barrier potential V 
are: 

POO 

E n (X) = (n+l/2)hu c + < n,X\V6(x)\n,X > = {n+l/2)hu c +V / dx p°JX,x) 

Jo 

(7) 

where p° n (X, x) is the electron density corresponding to the solution of eq.(3) 
without the barrier. It is given by 



p° n (X,x)=p^ n (x-X) 
1 



p L n an {x) = ' , - exp 



x 2 



H 2 n (x/l) (8) 



P 

where H n (x) are the Hermite polynomials of the order n. The integral can 
be performed to obtain analytic expressions for the first order corrections of 
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the energies: 

AE = V [1/2 + 1/2 Erf(X//)] 

X/l 



AE 1 = V 



1/2 + 1/2 Erf (X//) - ^= exp (-X 2 /f 

2a/7T v 



(9) 



for the first two levels. Differentiating eq.(7) with respect to X one obtains: 

vrfr(x). (io) 



dE n (X) . , 



dX 

The inverse density of states of the n th Landau branch is proportional to 
9E q^ ■ Now, since the electron density, eq.(8), has n nodes, perturbatively, 
there are precisely n infinitely flat regions of the dispersion relation. These 
will be called edge plateaus. For example, Ei(X) has such region around 
X = 0. Generally, the zeros of the Hermite polynomials thus determine the 
plateau centers. 

The gap between the ground and the first excited Landau branches does 
not remain Tioj c . For X > it gets reduced and reaches its minimum value 
hu c - ^j^V at X = l/y/2. 

These two features turn out to be of general nature rather than just an 
artifact of the perturbation theory. As we show next, they become even more 
pronounced for higher barriers. 

3.2 Spectrum. 

In order to solve the quantum mechanical problem eq.(3) with the rectan- 
gular potential for arbitrary height V beyond perturbation theory, the usual 
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quantum mechanical matching of the wave function has to be performed 
along the boundary line x — 0. This was done for the finite barrier in [?] (in 
which, however, the edge plateaus have not been noted). 

The wave function matching for (f> n> x(x) has to be performed at a single 
point x — 0. It is convenient to shift the origin of the coordinate system to X. 
In natural units of magnetic length the new variable x' is: x' = ^f{x — X). 
With energy expressed in units of Landau spacing hu c , v n = E nt x/{hu c ) — |, 
eq.(3) takes a form 

d 2 



_ _ V9{x > + y/2X/l) + (v n + I) 



<p n ,x(x') = (11) 



dx' 2 

which is the differential equation defining the parabolic cylinder functions 
[?, ?]. The two linearly independent solutions D Un (—x') and D Vn ^y(x') satisfy 
asymptotically the conditions of rapid decrease for x' — > — oo and x' — > +oo, 
respectively. 

The matching of the logarithmic derivatives at X' = —\/2X/l (or equiv- 
alently the condition of zero Wronskian in the expression for the Green's 
function of the system) gives: 

D' Vn „ v {x>)\ xl=x , D' Un (-x')\ x , =x , 

This determines the energy levels v n as functions of the position X. The 
equation (12) was solved numerically using a simplified form (obtained from 
the well known relations expressing derivatives of the parabolic cylinder func- 
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tions [?]): 

D Un+1 {X')D Un _ v {-X') + D Un _ v+1 {-X')D Un {X') = (13) 

The dispersion relations for electrons E n (X) (Landau branches) in the 
vicinity of the barrier have been evaluated numerically. Results for V = hu c 
and V = 5hu c are given in Fig.l and Fig.2[|. The ground Landau branch 
behaves as expected. It starts at its left asymptotic value v = and gradually 
rises to the right asymptotic value V + v. The spectrum of excited Landau 
branches above the top of the barrier is rather surprising, however. Instead of 
smooth transition from one asymptotic region (X << 0) to another (X » 
0), a steplike rise is obtained. Two unexpected features can be clearly seen: 

(i) Edge plateaus. 

The n th Landau branch E n (X) has n very pronounced edge plateaus. We 
call these edge plateaus in order to differentiate them from those in the 
asymptotic regions (the bulk plateaus). Density of states at the edge plateau 
energies is strongly enhanced. 

(ii) Landau gap reductions. 

The n th Landau branch almost touches the (n — l) th one n times. Unlike in 

the case of smooth confinement in which all the Landau branches follow the 

underlying potential (], the energy gap between two neighbouring branches 

exhibits regions of significant reduction, in particular for higher barriers. 
4 Results for a high barrier V = 5 huj c were reported earlier [?]. 

5 Similar behaviour is observed in the opposite extreme case of abrupt and infinite 
barrier [?]. 
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Following a Landau branch edge plateaus and the gap reduction regions 
alternate. First we would like to understand the origin of the edge plateau 
phenomenon. 

3.3 Why there are infinitely flat plateaus for the step 
potential? 

The plateaus in Fig 1 and 2 look, in fact, very flat. How flat? In order to 
see this, we calculate the density of states D(E). As it is well known [?] the 
density of states in ID system is just an inverse of the derivative 



DJE) = - 



dE n (k y ) 



2nP 



dE n (X) 
dX 



(14) 



Fig. 3 shows the density of states of the third excited n = 3 Landau level for 
V = 5hu c . The three peaks around v — 5.7, v — 6.4 and v — 7.1 indicate 
infinite density of states. In Fig. 4, full density of states from all the branches 
below Fermi level E F = 3.5 hu c is plotted. 

Now we will show that this is indeed the case for any rectangular potential. 
We apply the Hellmann - Feynman theorem to Hamiltonian eq.(ll), again 
with respect to X, to obtain: 

dE n (X) 



dX 



-< <p nX \V5(x' + X)\<j) nX >= V\<p nX (x' = -X)\' 



(15) 



Returning back to the variable x the result is 

dE n {X) 



dX 



V<j> 2 nX (x 



0) 



(16) 
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This is the square of the (real) wave function at the interface. Note that 
from eq.flTED we can obtain the monotonicity of i? n (X)[] 

THI FILE ARRIVED TRUNCATED. THE AUTHORS DIDN'T CARE. 
NEITHER DO WE. 



6 It is convenient to use eq.(|l6|) to normalize the wave functions. This circumvents the 
numerical integration. 
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